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ABSTRACT 

We investigate how observations of strong lensing can be used to infer cosmological parameters, in particular the equation of state 
of dark energy. We focus on the growth of the critical lines of lensing clusters with the source redshift as this behaviour depends on 
the distance-redshift relation and is therefore cosmologically sensitive. Purely analytical approaches are generally insufficient because 
they rely on axisymmetric mass distributions and thus cannot take irregular critical curves into account. We devise a numerical 
method based on the Metropolis-Hastings algorithm: an elliptical generalization of the NFW density profile is used to fit a lens model 
to an observed configuration of giant luminous arcs while simultaneously optimizing the geometry. A semi-analytic method, which 
derives geometric parameters from critical points, is discussed as a faster alternative. We test the approaches on mock observations 
of gravitational lensing by a numerically simulated cluster. We find that no constraints can be derived from observations of individual 
clusters if no knowledge of the underlying mass distribution is assumed. Uncertainties are improved if a fixed lens model is used for 
a purely geometrical optimization, but the choice of a parametric model may produce strong biases. 

Key words, gravitational lensing: strong - clusters: general - cosmological parameters - dark energy 



1. Introduction 

In the widely accepted 'concordance model' of cosmology a cos- 
mological constant accounts for more than 70 per cent of the 
overall energy density in the universe. It is considered an impor- 
tant feature mainly because it can account for both the spatial 
flatness and the accelerated expansion of the universe. Still there 
is no physical explanation that is generally considered satisfying, 
and the observational evidence can be reproduced in models that 
instead introduce dark energy, characterized by negative pres- 
sure. Observations of gravitational lensing may help determine 
its equation of state. Lensing phenomena are sensitive to the ge- 
ometry of the cosmological background since the appearance of 
an image depends on the distances between source, lens and ob- 
server. If we can obtain information about these distances from 
observations, we can relate them to redshifts to constrain space- 
time curvature, which is governed by cosmological parameters. 

In this paper we study strong gravitational lensing. Since a 
mass distribution needs to feature very high densities to act as 
a strong lens, galaxy clusters are suitable subjects of investiga- 
tion. The critical lines of such lenses grow with the redshift of 
the source. The methods discussed here aim to infer geometri- 
cal information from observations of giant arcs, which trace the 
critical lines, and consequently constrain the equation of state 
of dark energy. An exploratory study of the concept was pre- 
sented by Meneghetti et al. (2005 ). Other authors followed simi- 
lar approaches to analyse individual clusters, e . g. Sereno (2 002) , 
Soucail et ah] ( |2004| >, |Gilmore & Natarajan| ( |2009| >. Moreover 
weak lensing can be studied with the same goal, as presented by 
Medezinski et aT] ( |201 1| , for instance. 

We give a short summary of the underlying theory in Sect. [2] 
In Sect. [3] we present analytic studies of cluster lensing. While 
they are restricted to axisymmetric lenses, they can provide us 
with estimates of the influence of the dark energy equation of 
state on strong lensing features. In Sect. [4] we turn to numeri- 



cal approaches and develop Markov chain Monte Carlo methods 
that aim to fit a set of parameters, characterizing the lens model 
and the geometry, to an observed image configuration. In Sect. 
[5] an alternative approach is presented, which infers the geom- 
etry from the observed scaling of the critical lines at different 
redshifts. Finally in Sect. [6] we discuss the performance of the 
various approaches and point out possible sources of errors. 



2. Theoretical background 

2.1. Cosmological model 

We assume that the universe is spatially flat and characterized by 
the Friedmann-Lemaitre-Robertson-Walker metric 

ds 2 = -c 2 df 2 + a 2 (f) [dw 2 + w 2 (d6 2 + sin 2 6d<f> 2 )] , (1) 

where t denotes the coordinate time, w is the comoving radial 
coordinate and Q and cp are the azimuthal and polar angles. The 
evolution of the scale factor a( f) is governed by the Friedmann 
equations, 

U = — p (2) 



and 



AnG 



3p 



3 \ P+ c> 



(3) 



We have omitted the curvature terms as well as those involv- 
ing a cosmological constant. Instead we consider a dark energy 
component with the equation of state p x = wp x c 2 , such that its 
density is given by 

p x (a) = p xfi exp ^-3 1 + W ^ a - da' j . (4) 
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To drive the expansion of the universe w < -1/3 is required. 
From Eq. (|2]) the Hubble function H = a/a can be derived: 



H(z) = H A Q m , (l + z) 3 + (1 - Q m ,„) exp 3 



1 + w(z') 



dz'i 77 



jo 1 + z' 

(5) 

Here the matter density parameter Q. m .o has been introduced. 
Since curvature is assumed to vanish and the radiation density 
is negligible in the epochs probed by gravitational lensing, the 
dark energy density parameter is Q Xj o = 1 - Q m ,o- 

2.2. Gravitational lensing 

We summarise only the key aspects of gravitational lensing 
here. A comprehensive review of the theory was provided by 
|Bartelmann| ( T2010l l. 

The relation between the angular positions of an image, 0, 
and its source, /?, is given by the lens equation: 

P = 6 - oe(0). (6) 

a is called the reduced deflection angle and can be written as 



(7) 
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Fig. 1. Growth of the Einstein radius # E of a singular isothermal 
sphere (cr = 1000 km s _1 ) with the source redshift z s for different 
values of the equation of state parameter w. In the lower panel 
the deviations from ACDM are shown. 



where the deflection angle a depends only on the lensing mass 
distribution. D ds is the distance from the deflector to the source 
and D s the distance from the observer to the source. Both are 
angular diameter distances, with the distance to redshift z given 
by the integral 



1+zJo H(z') 



(8) 



Reduced deflection angles thus depend on the lens and source 
redshifts. 

In the thin-screen approximation the lensing potential \f/ is 
defined as the Newtonian potential projected on to the lens 
plane, such that its gradient is the reduced deflection angle, 
a{6) = Vgi//(0). Image distortions are quantified by the Jacobian 
J of the lens mapping (7y = d/3i/d6j), specifically by the conver- 
gence 



K{0) 



and the shear 
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(9) 
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The convergence is proportional to the Laplacian of the lens- 
ing potential, V 2 i/r(0) = 2k{0). It measures the projected surface 
mass density £ of the lens in units of the critical density: 



1(0) 



4nG D d D ds 



(11) 



Critical lines are formed by those points in the lens plane where 
the lens mapping is singular, det J = (1 - k) 2 - y 1 — = 0. As 
the lensing potential changes with the source redshift, so does 
the shape of the critical curves. 



3. Analytical approaches 

3. 1 . The singular isothermal sphere 

We first examine the lensing behaviour of the singular isothermal 
sphere profile. In this model the density of a galaxy cluster is 
described by the function 



P(r) = 



cr 2 
2nGr 2 ' 



(12) 



where cr is the velocity dispersion of the cluster members. Such 
a lens is characterised by the Einstein radius 



ft, = 4;r 



(cr\ 2 Dd, 
\c) A 



(13) 



Figure [T] shows the dependence of the Einstein radius on the 
source redshift in different dark energy cosmologies for a cluster 
with velocity dispersion cr = 1000km s~' at redshift Zd = 0.3. 
In addition to the ACDM scenario w = -1, we choose to ex- 
amine the behaviour for a constant equation of state parameter 
w = -0.7 and an extreme case of dark energy, a phantom model 
with w = -1.3. While the Einstein radius itself takes values of up 
to 25" for the redshift range considered, differences between the 
cosmologies are less than 1" and largest if the source is located 
close to the cluster. 

To obtain the angular diameter distance ratio D^/D^ from a 
single measured Einstein radius 8e the velocity dispersion cr of 
the lens must be known. We can eliminate it from our analysis if 
we study the growth of the critical curve with the source redshift, 
comparing two Einstein radii #ei, &e2 at different redshifts z s \, 
z%2- The resulting ratio 



j. _ 8e2 _ D&s2 D s \ 

c? E i D S 2 Ajsl 



(14) 



appears frequently in our analyses, and we label it the geom- 
etry factor. Via this factor, cosmological parameters determine 
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Fig. 2. Growth of the geometry factor / with the source redshift 
z S 2 for two different dark energy cosmologies. The fixed refer- 
ence redshift is z s i = 1.0 or z s i = 0.7 respectively. 

lensing properties of an object. It should be noted that the pro- 
portionality between the geometry factor and the Einstein radius 
is peculiar to the SIS model. For a fixed redshift z s i, the geom- 
etry factor increases with the redshift z S 2 of the second source, 
rising steeply behind the lens and flattening out for high red- 
shifts. In Fig.|2]this behaviour is displayed for a lens at redshift 
Zd = 0.3 and fixed sources at z s i = 1.0 or z s i = 0.7 respec- 
tively. Differences between the two cosmologies shown are most 
prominent if the two sources are at a high distance from each 
other, since this configuration corresponds to a large 'lever arm'. 
For a cosmological analysis it is therefore the most convenient 
to study pairs of arcs in which one source lies closely behind the 
lens and the other at a much higher redshift. 

3.2. The NFW profile 
The NFW profile 



P(r) 



(15) 



r(r s + r) 2 

is arguably the most commonly used parametric model for the 
matter distribution of dark matter haloes. For its Einstein radius 
no analytic expression exists, but Bartelmann (1996) provided 
the convergence: 



k(x) = 



2k s 



1 



1 - 
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- arctanh 
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(16) 



x — r/r s is the dimensionless radial coordinate in the lens plane 
and k s = p s r s S^. ; note that Eq. ([16) is valid only for x < 1, i. e. 
inside the scale radius, which encompasses the strong lensing re- 
gion. We start our analytic approach from a property of axisym- 
metric lenses, namely the condition that the mean convergence 
inside the tangential critical line of radius 0e has to equal unity: 



i r fe 

nOl Jo 



2n9d0K(6). 



(17) 



If we set up this equation for two different source redshifts z s i 
and z S 2, rearrange and divide the two equations, we arrive at the 
form 

Aisi D s2 _ #ei 8 (*E2) 



Ai £> d s2 6% 2 8(xm) 
The function g(x) represents the integral 



(18) 
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(19) 



solved by Bartelmann (1996). Again we have used dimension- 
less coordinates in the lens plane, i. e. x El2 = D^6 El2 /r s . The 
left-hand side of Eq. ( fTS) is a function of the source redshifts, 
while the right-hand side depends on the Einstein radii for those 
redshifts; both sides also depend on cosmological parameters 
since these determine the distance-redshift relation. In this form 
the equation admits a simple graphic solution: plotting each side 
against the equation-of-state parameter w, the intersection of 
both curves marks the true value. 

To test whether the relation can be exploited that way, we 
consider as an example a halo of mass M = 1.0 x 10 15 hr l Mq 
and scale radius r s = 310/^' kpc at a redshift of za = 0.3 and 
sources at redshifts z s i = 0.7 and z S 2 = 3.0. The Einstein radii for 
those redshifts are E1 = 7'.' 9 and E2 = 21'.' 9 in ACDM cosmol- 
ogy. In reality, we have to resort to estimates of the scale radius, 
such as best-fitting values, which introduce errors. Einstein radii 
have to be determined from the positions of observed arcs. While 
the latter can be expected to trace the critical lines of the clus- 
ter, uncertainties in the deduced Einstein radii may be at least as 
large as their widths and thus of the order of 1". 

Figure [3] shows the influence of such errors on Eq. (18) , In 
the first plot, the correct values for the Einstein radii are used, 
but the estimate for the scale radius is too high or too low re- 
spectively, while in the second plot one of the Einstein radii is 
underestimated. In each calculation, the correct values are used 
for all but the specified quantity. The left-hand side of the equa- 
tion is independent of such errors since, as stressed before, it 
depends only on the source redshifts, of which exact knowledge 
can safely be assumed here. 

Inaccuracies in the scale radius estimate do not have a very 
large effect on the result for the equation of state parameter, shift- 
ing the value of w by Aw ~ 0.1 in this example. The determina- 
tion of the critical line poses a larger problem, with an error of 
A0 E ~ 1" in the Einstein radius translating into a deviation of 
Aw ~ 0.6. Another caveat is given by the fact that galaxy clus- 
ters often possess significant ellipticity. Based on our analysis 
so far it seems unlikely that these simple, spherically symmetric 
profiles provide a sufficient approximation, if only for the fact 
that it is unclear how a robust measure of the Einstein radius can 
be obtained from arc positions in the case of non-circular critical 
lines. 

4. Markov chain Monte Carlo methods 

4.1. The lens model 

To avoid the assumption of axisymmetry, we focus on an ellipti- 
cal generalization of the NFW profile used by Come rford et al.| 
(2006). This model is based on six parameters that character- 
ize the cluster constituting the gravitational lens: the coordinates 
(x c j c ) of its centre, the scale convergence k s , the scale radius 
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Fig. 3. Left-hand side (LHS) and right-hand side (RHS) of Eq. 
( p~8] > plotted against the equation of state parameter w for various 
estimates of the scale radius r s (upper panel) and Einstein radii 
#ei, #E2 (lower panel). For each curve only the specified quan- 
tity is changed from its true value. The thick line in each panel 
corresponds to the correct values r s = 310 h kpc, 9%\ = 7'.' 9, 
fta = 21'.'9. 



r s , the ellipticity e and the position angle (f>. To calculate the de- 
flection angle field the coordinate frame is shifted and rotated in 
such a way that the cluster centre determines the origin and the 
coordinate axes coincide with the axes of the ellipse. Then an 
elliptical radius is introduced, 



p = -v/x 2 (l -€) + 



r 



1 -e 



(20) 



To obtain the deflection field, the NFW lensing potential if/ is 
evaluated at this radius and a is calculated by differentiation 
with respect to the coordinates in the original frame. The conver- 
gence k and shear y can be computed by further differentiation. 
Deflection angles obtained in this way are valid only for a refer- 
ence source redshift z r . For each additional source i at a different 
redshift z s i, values have to be multiplied by the geometry factor 



fi 



£>Afa,Zsi) £>A(0,Zr) 

D A (0,z si ) D A (z d ,z r )' 



(21) 



Given an image point originating from source i, the corre- 
sponding source point ft is located using the lens equation 



(22) 



with the appropriately rescaled deflection angle fia entering. By 
scanning a grid on the lens plane for points that fulfil the lens 
equation for this source position, all image points are located. 



4.2. The chi-square 

To quantify how well a set of parameters describes the observed 
lensing effects, we follow Comerford et al. (2006) in introducing 



a ^-function that has three contributions, x 



-■x\+xl+x\- 



We demand that the first should measure the extent to which 
the observed images can be reproduced. Let N be the number 
of data points with coordinates (xi,yi) and (mj, vy) the coordi- 
nates of the predicted image points. For each data point, the 
closest image point {u c \^, v c y) is identified, leading to 



-y 

N 4^ 



{Xi - M c l,i) + CVi - y cl,i) 



(23) 



This form assumes that data points are distributed around the 
predicted image points in a Gaussian fashion, with a standard 
deviation of cr ; . 

- The second contribution should test whether the predicted 
images match the data. For each predicted image point (out 
of M overall), the closest data point is identified. The chi- 
square contribution is then given by 



xl 



-7 

1=1 



(Mi - Xcu) + (vi - y cU ) 



(24) 



If the model gives rise to any additional image points far 
from the data, it affects this term. Obviously this contribution 
has the same form as the first, but the roles of data and image 
points are reversed. 

Finally, the size of the sources is taken into account. For each 
of the A^ s sources, the centre {p,,qi) is determined by com- 
puting the mean of either coordinate. Then the mean squared 
distance from the centre is calculated, averaging over the P; 
points assigned to the j'-th source. For a tolerated source size 
cr s we take 



xl 



1 



z 

1=1 



i A {pur_ 

7=1 



Pi)' + [qui - q~if 



(25) 



as the final chi-square contribution. Again it is assumed 
that the distribution of points belonging to the same source 
around its centre is Gaussian. This last contribution is in fact 
crucial because neglecting it would mean that source shapes 
and sizes are arbitrary - in that case, any image configuration 
could easily be reproduced with a lens of zero mass and the 
points placed in the source plane perfectly matching the dis- 
tribution in the image plane. Requiring instead that sources 
are small and compact is therefore a strong constraint. 

We find that a reasonable choice is cr ; = l'/O and <x s = 075, 
emphasizing the source size. Note that the behaviour of the al- 
gorithm is determined by the ratio between the two parameters, 
not their absolute sizes. 



4.3. The Metropolis-Hastings algorithm 

We aim to fit a set of lens parameters and geometry factors to ob- 
served arcs. Because of the large number of free parameters (six 
for the lens model and one for each pair of sources) we avoid us- 
ing a simplex algorithm to minimize the chi-square. Instead we 
generate Markov chains according to the Metropolis-Hastings 
algorithm. Given a set of parameter values x in \ a random point 
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Jn+l) 



is picked from the parameter space and accepted with the 



transition probability 



t(x M 



) = 



p{x^)q{x in) x (n+1) ) 
l ' p(x^) q(x ( "+V -> xf- n >) 



(26) 



p(x) is the target probability distribution. We choose the likeli- 
hood L = Ne\p(~x 2 /2) (where N accounts for proper normal- 
ization), so that the sample is concentrated on regions of high 
likelihood. By means of the proposal density q (x — > y) the step 
sizes can be limited or parameter ranges set. 

4.4. The Skylens simulator 

To test the algorithm on mock observations, we use SkyLens, a 
ray-tracing code pre sented by |Meneghetti et al.|(2008|l. It was 
used for instance by |Merten et al. ( 2009 1 and Meneghetti et al. 
(2010). In short, the program generates a distribution of back- 
ground galaxies, based on a set of real galaxies decomposed into 
shapelets. Positions in a specified field of view and orientations 
are randomly selected. If desired, all sources can be placed at a 
fixed redshift: this feature is very useful for studies like ours that 
are based on observing the change of lensing properties with the 
source redshift. Observational effects are added to the lensed im- 
age, including the sky background, photon noise and seeing as 
well as instrument noise. The lensed images can be convolved 
with point spread functions, which are available for several tele- 
scopes. 

We choose to study a numerically simulated cluster la- 
belled gl, taken from a sample of hydrodynamical simula- 
tions by Saro et al. (2006i. It was obtained from a dark mat- 
ter simulation by Yoshida et al. (2001 1 and re-simulated with 
added baryonic effects at a higher mass and spatial resolution 
using Gadget-2 (Springel 2005). The cluster has a mass of 
M200 = 1.14 x ltP h~ l Mq and a best-fitting scale radius of 
r s = 0.310/i -1 Mpc. Principal axis ratios are b/a = 0.64 and 
c/a = 0.57 and the orientation of the main axis relative to the 
coordinate axes of the simulation box is given by the angles 
X = 33.3°, G y = 57.4° and G z = 96.1°. The cosmological pa- 
rameters used in the simulation are Q A o = 1 - O,,, = 0.7 
(with Q/, o = 0.04) and h = 0.7. Detailed explanations can be 
found in papers about other studies using these simulations (e. 
g. |Dolag et alT] ( |2005| ); |Puchwein et al.| ( |2005| )). D eflection an- 
gle maps were computed for different projections bylMeneghetti 
|et al.| ( |2008"] ), placing the cluster at redshift z d = 0.2975. 

We created mock observations of the mass distribution pro- 
jected along the z-axis for the Advanced Camera for Surveys 
(ACS) on HST and obtained 9 giant luminous arcs at redshifts 
Z = 0.7, 1.0,2.0,4.0. 

4.5. Results 

We choose flat prior distributions for the lens parameters, confin- 
ing parameter values to fixed intervals. If no assumptions at all 
are made about the lens, it is difficult to obtain information about 
the geometry due to the degeneracies involved; e. g. trends in 
the scale convergence and geometry factor can compensate each 
other to some extent. Moreover, confining the parameters makes 
the parameter space smaller and accelerates its exploration. For 
the geometry factors we do not explicitly exclude any region 
from the beginning and limit only the step size as the algorithm 
is designed such that the chain should generally move 'in the 
right direction' regardless of the starting point. In practice, start- 
ing values for the geometry factors must still not be set too high 



unless care is taken to ensure the correct calculation of the tiny 
likelihoods and their ratios in particular. Generally starting val- 
ues between and 5 and search radii between 0.2 and 0.5 lead 
to reasonable burn-in phases and acceptance rates. 

We consider only one pair of sources at a time: one located 
at the reference redshift for the lens parameters, the geometry 
factor of which consequently has the value 1 .0 and is not varied; 
and a second source at a different redshift. Including the refer- 
ence redshift helps break the degeneracy between the geometry 
factor / and the scale convergence k s . Otherwise changes in ei- 
ther quantity can be absorbed in the other, since deflection angles 
are proportional to the product f ■ k s . We take the lowest source 
redshift of z s = 0.7 as the reference redshift and follow the pro- 
cedure described above to produce Markov chains for each ge- 
ometry factor. From these samples we compute likelihood dis- 
tributions, marginalizing over the lens parameters. Histograms 
for the distributions are presented in Fig. |4] In each plot, the 
true value of the geometry factor, i. e. the value in the ACDM 
cosmology assumed in the simulation of the images, is marked. 
The locations of the likelihood peaks generally agree quite well 
with the true values of the geometry factors. Yet as pointed out 
in Section [3] geometry factors vary very little between different 
dark energy cosmologies (cf. Fig.[2|. As w ranges from the rather 
extreme scenario w — -2 to w = -0.3 for instance, geometry 
factors vary by less than 0. 1 or roughly 6 per cent (depending on 
the redshift). On the other hand, the widths of the distributions, 
which for simplicity we quantify using the standard deviation of 
the best-fitting Gaussian (despite the skewness), are roughly 8 
per cent of their respective mean. The problem lies in the fact 
that the lens model can easily be adjusted to react to any small 
change in the geometry. To infer cosmological information, how- 
ever, a much better 'resolution' is needed. 

To address this problem we explore how the method behaves 
if the lens parameters are kept constant. Ideally, no geometrical 
assumptions should be made in choosing the lens parameters. In 
the absence of independent information from effects other than 
strong lensing, this means that only images at one redshift may 
be used to fit a model. We attempt this using the code presented 
by Comerford et al. ( |2006| ). However, we find that it is not possi- 
ble to obtain a reliable parameter set in this way as fit results vary 
very strongly with different choices of starting values. As we are 
nonetheless interested in the performance of the algorithm for a 
fixed lens model, we perform a fit including arcs at three differ- 
ent redshifts. Again we study the likelihood distributions for the 
geometry factors, shown in Fig. [5] Compared to the full variation 
of both lens parameters and geometry factors, the distributions 
do appear considerably more narrow, with widths of 2-3 per cent 
of the mean. However, deviations from the ACDM values are 
still not satisfactory. In most cases, the likelihood peaks are lo- 
cated at geometry factor values that are higher than in ACDM. 
Comparing the results for several arcs at the same redshift, we 
note that the likelihood distributions do not seem consistent in 
that they do not appear to favour the same geometry factors. 

In the following we investigate whether our choice of the 
lens model can account for such deviations. 



4.6. Influence of the lens model 

To check whether the strong biases observed in the likelihood 
distributions presented above are indeed caused by insufficient 
knowledge of the lens model, we repeat the procedure using 
arcs produced by an analytic deflection angle map rather than the 
simulated cluster. In order to be able to carry out calculations an- 
alytically, we define the deflection angle field by a simple power 
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The dashed vertical line in each plot indicates the value of the each plot indicates the value of the geometry factor in a ACDM 



geometry factor in a ACDM cosmology. 



cosmology. 



law: 



a(x) — aox p 



(27) 



x — r/ro is the dimensionless radial coordinate in the lens plane; 
a is also dimensionless and denotes the deflection angle relative 
to the angle set by the reference scale ro. The latter is arbitrary 
but has to comply with ao = a(x = 1). The map is therefore 
characterized by two parameters. 

We choose a scale of ro = 0.3/r'Mpc and set ao = 0.26 
and p = 0.36, such that the deflection angle field approximately 
follows that of an NFW halo of the same scale radius ro and the 
scale convergence k s = 0.2; the behaviour for both the power law 
and the specified NFW profile is shown in Fig. [6] We demand 
that the field should describe the deflection angles for a reference 
source redshift of z s i = 0.7 with the lens located at za = 0.3. 
In SkyLens simulations, this configuration produces two giant 
luminous arcs originating from the same source at redshift z S 2 
4.0. 

The MCMC method again provides likelihood distributions 
for the geometry factor. Since we want to study the influence 
of the assumed lens model, we run it several times, varying the 
parameters ao and p. ro is kept fixed; changing its value has the 
same effect as changing ao. As seen in Fig. [7] the modifications 
clearly move the likelihood peaks, but they hardly change the 
shape of the distributions. Based on our exact knowledge of the 
underlying mass profile in this case we can attempt theoretical 
predictions: The deflection field described by Eq. |27) leads to a 
convergence of 



for the reference redshift; for other redshifts this has to be 
rescaled by the geometry factor /. Now we resort to Eq. fP7] i 
again to compute the critical line. The mean convergence inside 
the radius x c of the critical line is 



(29) 



(30) 



1 = faoxl , 
which leads to an Einstein radius of 

i 

x c = (/ao) ■ 

For our choice of parameters, this gives a radius of x c i = 0.122, 
corresponding to 11'.' 8, for the reference redshift z s i - 0.7. The 
geometry factor for redshift z S 2 — 4 is / = 1 .5925, consequently 
the Einstein radius should be x C 2 = 0.252 (24'.'4), in good agree- 
ment with the observed arcs. 

Conversely, given a profile (ao,/?) and an Einstein radius x c 
we can calculate a geometry factor /, solving d29) accordingly: 



/--!■« 

a 



(31) 



k(x) — ^— — aox p 1 



(28) 



This enables us to test how a wrong choice of parameters affects 
the deduced geometry factor and compare the results to the like- 
lihood peaks. In the histograms in Fig. [7] the peak locations ex- 
pected from the analytical estimation are marked along with the 
correct values. They agree remarkably well with the actual like- 
lihood maxima. It seems that the shift in the distributions away 
from the true geometry factor can indeed be accounted for by an 
erroneous choice of mass profile parameters. Since such a choice 
does not change the features of the distributions, it should gen- 
erally be difficult to distinguish it from other cases with better 
choices and remove the effect. 
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Fig. 6. The dimensionless deflection angle a as a function 
of the radius r in the lens plane, for a power law a(r) = 
0.26[r/(0.3/z _1 Mpc)] 036 and an NFW halo with scale conver- 
gence k s = 0.2 and scale radius r s = 0.3 h~ l Mpc. 
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Fig. 7. Likelihood distributions L for the geometry factor / of 
an arc produced at redshift z.a — 4 by a power law deflection 
field. Coefficient ao and slope p are fixed at the indicated values; 
true values are a = 0.26, p = 0.36. The thick dashed line marks 
the true geometry factor, thinner lines indicate predictions from 
analytic calculations for the respective parameter choice. 



5. Semi-analytic calculation of geometry factors 

We test a final approach that is based on the same parametric 
profile but aims to reproduce critical points rather than images. 
While it seems that arc configurations can often be predicted by 
a variety of lens models in different geometries, to some extent 
owing to the freedom in the source shapes and the position of 
images relative to the critical line, the shape of the critical curves 
itself may be more demanding to reproduce. 



Table 1. Geometry factors / computed for three source redshifts 
z s from the chi-square minimisation, with a reference source red- 
shift of either z r = 0.7 or z r = 1.0, and the theoretical ACDM 
values. For a reference source redshift of z r = 1.0, the geometry 
factor for source redshift z s = 1 -0 is / = 1 by definition. 



Zr 


Zs 


/ (comp.) 


/ (ACDM) 




1.0 


1.2334 


1.2168 


0.7 


2.0 


1.5256 


1.4634 




4.0 


1.6058 


1.5839 


1.0 


2.0 


1.2369 


1.2027 




4.0 


1.3019 


1.3017 



We consider a lens model with convergence kq and shear 

y° ~ sl^o l + y\ 2 f° r a source redshift z r . For a second source 
redshift z s characterized by the geometry factor / the con- 
vergence at any point (x,y) in the lens plane has the value 
k(x, y) = /kq(x, y); the shear also scales with the geometry factor, 
y(x,y) = fyo(x,y). On the critical lines 



detJ(x,y) = [1 - K(x,y)Y - y 2 (x,y) = 



(32) 



has to hold. Given a set of N points {Xi,yi) known to lie on the 
critical line at redshift z s , we therefore take the function 



X 



1 N 2 

jj J] {[1 - AofeyO] 2 - fVoix-.ydf (33) 



to measure the agreement between the model and the data for 
the geometry factor / and demand that its minimum should de- 
termine the best-fitting geometry factor. To apply this concept 
we first have to obtain sets of critical points. Ideally an estimator 
point is obtained by marking the brightness saddle point in an 
arc that appears to be formed by two merging images. In arcs 
exhibiting no such structure the brightest point can be chosen 
instead. It is useful to add more points nearby to mark the pre- 
sumed direction of the critical line in this point since this makes 
for more stringent constraints. Multiple image systems can also 
be taken into account. If several arcs are observed at the same 
redshift at different positions in the lens plane, this translates 
into more information on the critical line and can help to con- 
strain ellipticity in particular. 

We test the method using our mock observations of the clus- 
ter gl. Since very few and only faint arcs are found at redshift 
z = 0.7, we include only redshifts z - 1, z — 2 and z = 4. Note 
that in principle we could fix the lens model for one of these 
redshifts by setting / = 1 in Eq. ( p3j ) and minimising with re- 
spect to the lens parameters. However it is all but impossible to 
determine a unique best-fitting model since a large number of 
combinations of lens parameters produces critical curves pass- 
ing through our estimators, some with scale radii several times 
as large as the value suggested by other mass profile fits for this 
cluster. To avoid this obstacle, we rely on the same set of fixed 
lens parameters as for the MCMC method. We then use estima- 
tor sets for all three redshifts to compute their geometry factors. 
Table[T]shows the values obtained in this way. 

Two different reference source redshifts are listed for the fol- 
lowing reason: The value for the scale convergence k s is origi- 
nally defined here for a source redshift of z r = 0.7. If we are con- 
fident that the model is valid for this redshift, we can consider 
the geometry factors with respect to the same reference redshift. 
If we do not trust the model however, we can restrict our analysis 
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Table 2. Influence of the scale radius r s and ellipticity e estimates on the geometry factor / in the semi-analytical method. In the 
third column, the previous parameter choice and results are given. Results in the remaining columns were computed with a changed 
value of either r s or e as indicated in the header. 



Zr 


Zs 


r s = 0.308 /j-'Mpc, e = 0.37 


r s = 0.317/j-'Mpc 


r s = 0.294 h-' Mpc 


e = 0.35 


e = 0.40 




1.0 


1.2334 


1.2150 


1.2631 


1.2707 


1.1782 


0.7 


2.0 


1.5256 


1.4991 


1.5685 


1.5800 


1.4452 




4.0 


1.6058 


1.5771 


1.6522 


1.6573 


1.5297 


1.0 


2.0 


1.2369 


1.2338 


1.2418 


1.2434 


1.2266 




4.0 


1.3019 


1.2980 


1.3081 


1.3042 


1.2983 



to the redshifts that actually appear in observations. Since 



/(z s = ZuZv = Zl) 



f(Zs 



Zl,Zt 



zo) 



f(Zs 



Z2,Z r 



zo) 



(34) 



for arbitrary redshifts zo,i,2 according to the definition of the ge- 
ometry factor (cf. Eq. (fT4]>), we consider ratios between the com- 
puted geometry factors. While for instance /(z s = 2.0, z r = 0.7) 
is the factor by which the input scale convergence has to be mul- 
tiplied for the critical line to match our estimator set for red- 
shift z s = 2, / (z s = 2.0, Zt = 1 -0) instead can be taken to quantify 
the scaling between our two estimator sets. In other words, the 
convergence is rescaled to reproduce the critical line at redshift 
z s = 1.0 and geometry factors for higher redshifts are considered 
with respect to that redshift. This has the advantage of eliminat- 
ing the influence of the scale convergence from our calculations. 
Naturally, the remaining lens parameters can still act as sources 
of errors. 

Table [Tj confirms these considerations to some degree: for 
the reference source redshift of the lens parameters, z r = 0.7, 
deviations of the geometry factors are roughly 4 per cent for a 
source redshift of z$ = 2 and 1 per cent for z s = 4, whereas 
considering the ratios they reduce to 3 per cent and less than 
0. 1 per cent respectively. The excellent agreement in the latter 
case is certainly to some extent coincidental. Taking our lack of 
knowledge of the precise mass model into account, uncertainties 
should be considerably larger than the deviation itself. 

To find out how our parameter choice affects the results, we 
simply rerun the program for different lens models. Having just 
described how to eliminate the scale convergence, we investigate 
the role of the scale radius r s and the ellipticity e. Table[2]lists the 
geometry factors computed for various cases in which changed 
values for either the scale radius or the ellipticity were used. 

For a reference source redshift of z r - 0.7, the geometry fac- 
tor results vary over a range corresponding to about 8 per cent 
across our examples of lens models. Changing the reference red- 
shift - that is, again considering only ratios - the scatter reduces 
to 1 per cent. In that case, however, all results calculated for 
z s = 2 overestimate the geometry factor without exception, due 
to errors in either the assumed mass model or our choice of crit- 
ical line estimators. It is also worth noting that in this example 
we used multiple arcs at the same redshift to derive estimators. 
In practice, fewer and fainter arcs, albeit at more redshifts, might 
be available, complicating the choice of reliable estimators. 

6. Conclusions 

We have investigated several approaches to cluster strong lens- 
ing to test its use as a probe of spacetime geometry, in particu- 
lar the dark energy equation of state. Generally this has proven 
challenging as changes in the observables induced by varia- 
tions in the equation of state are very small. The methods that 
we have explored make use of the growth of the critical lines 



with the source redshift. We have tried to constrain a cosmo- 
logically sensitive ratio of angular diameter distances, using first 
a Metropolis -Hastings algorithm constructed to fit an observed 
image configuration and secondly a semi-analytic approach that 
optimizes the distance ratio such that a fixed lens model repro- 
duces a given set of critical points. To test our methods we have 
studied strong lensing by a numerically simulated cluster, for 
which we have created mock observations. 

Modelling uncertainties have turned out to be the most im- 
portant source of errors in our study. The simplest models for the 
lensing mass distributions, such as the SIS and the NFW profile, 
assume spherical symmetry. They offer the advantage of admit- 
ting a largely analytic analysis and thus help us estimate geo- 
metric effects, but they are generally not well suited to a study of 
lensing by a real galaxy cluster as in most cases deviations from 
the assumed symmetry are too large. 

Models that include ellipticity provide a more accurate de- 
scription of the mass distributions. The larger number of param- 
eters, however, produces degeneracies. If all model parameters 
are treated as free, any cosmological sensitivity is masked since 
the background geometry and profile properties can influence the 
lensing behaviour in a similar way. An attempt to optimise the 
lens along with the geometry for an individual cluster does not 
result in meaningful constraints on the distance ratio, but instead 
admits a wide range of realisations of dark energy. Constraints 
can be narrowed down somewhat if the lens model is fixed and 
only the geometry is optimised, but care must be taken in the 
choice of the model. While the influence of the overall mass as- 
sumed is weak if we choose to consider only the scaling of the 
critical curve between two redshifts, the remaining parameters 
can still create biases. Ideally, no geometric assumptions should 
be made to fix the lens model, yet information from arcs at a 
single redshift alone is not sufficient to constrain it. Independent 
information from other observations could be used to address 
the problem, but it is generally difficult to obtain for the clus- 
ter core. Meneghett i et al.| ( [20 05) suggested that the position of 
the brightest cluster galaxy (when present) or the centre of X- 
ray emission could be referred to to determine the cluster centre, 
but they pointed out that errors of several arcseconds can occur. 
Additional constraints could be derived from the galaxy veloc- 
ity dispersion, weak gravitational lensing or X-ray temperature 
profiles. 

It is also worth investigating whether the restriction to ana- 
lytical profiles can contribute significant errors. Galaxy clusters 
are generally 'lumpy', which gives rise to the question if an ellip- 
soidal NFW profile nonetheless describes the lensing properties 
well enough to permit cosmological conclusions. If even the best 
fitting NFW profile is insufficient, biases are to be expected. An 
interesting approach to cluster strong lensing that forgoes the 
use of the NFW profile was presented by Zitrin et al. (2009). 
Assuming that mass follows light in a cluster, they assigned a 
power law profile to each visible cluster galaxy, scaled by the ob- 



8 



B. Zieser and M. Bartelmann: Constraints on Cosmological Parameters from Strong Gravitational Lensing by Galaxy Clusters 



served brightness, and smoothed the resulting distribution. Note 
that they based their method on multiple images. While we have 
not considered this effect in our study, it provides additional con- 
straints and should therefore be a useful inclusion. 



References 

Bartelmann, M. 1996, A&A, 313, 697 

Bartelmann, M. 2010, Classical and Quantum Gravity, 27, 233001 
Comerford, J. M., Meneghetti, M., Bartelmann, M., & Schirmer, M. 2006, ApJ, 
642, 39 

Dolag, K., Vazza, F., Brunetti, G., & Tormen, G. 2005, MNRAS, 364, 753 
Gilmore, J. & Natarajan, P. 2009, MNRAS, 396, 354 

Medezinski, E., Broadhurst, T., Umetsu, K., Bem'tez, N., & Taylor, A. 2011, 

MNRAS, 414, 1840 
Meneghetti, M., Jain, B., Bartelmann, M., & Dolag, K. 2005, MNRAS, 362, 

1301 

Meneghetti, M., Melchior, P., Grazian, A., et al. 2008, A&A, 482, 403 
Meneghetti, M., Rasia, E., Merten, J., et al. 2010, A&A, 514, A93 
Merten, J., Cacciato, M., Meneghetti, M., Mignone, C., & Bartelmann, M. 2009, 
A&A, 500, 681 

Puchwein, E., Bartelmann, M., Dolag, K., & Meneghetti, M. 2005, A&A, 442, 
405 

Saro, A., Borgani, S., Tornatore, L., et al. 2006, MNRAS, 373, 397 
Sereno, M. 2002, A&A, 393, 757 

Soucail, G., Kneib, J.-P, & Golse, G. 2004, A&A, 417, L33 
Springel, V. 2005, MNRAS, 364, 1105 

Yoshida, N., Sheth, R. K., & Diaferio, A. 2001, MNRAS, 328, 669 
Zitrin, A., Broadhurst, T., Umetsu, K., et al. 2009, MNRAS, 396, 1985 



9 



